Solutions to the Boltzmann equation 
in the Boussinesq regime. 
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Abstract We consider a gas in a horizontal slab, in which the top and bottom walls 
are kept at different temperatures. The system is described by the Boltzmann equation 
(BE) with Maxwellian boundary conditions specifying the wall temperatures. We study the 
behavior of the system when the Knudsen number e is small and the temperature difference 
between the walls as well as the velocity field is of order e, while the gravitational force is 
of order e^. We prove that there exists a solution to the BE for t G (0, f) which is near a 
global Maxwellian, and whose moments are close, up to order to the density, velocity 
and temperature obtained from the smooth solution of the Oberbeck-Boussinesq equations 
assumed to exist for t <i. 

1. Introduction. 

In the study of thermal convection phenomena the following system plays a paradigmatic 
role: a viscous heat conducting fluid between flat horizontal plates with the lower plate 
maintained at a temperature greater than the upper one. When the temperature difference 
between the plates is small, the stationary state is one in which the fluid is at rest with a lin- 
ear temperature proflle. When the temperature difference is made larger, the gravitational 
buoyancy force acting on the light, higher temperature fluid below, overcomes the effects 
of viscosity and a new stationary state of thermal Rayleigh-Benard convection sets in. In 
typical experimental situations the variations of temperature and density are small and the 
system is described in the Boussinesq approximation under which the Navier- Stokes equa- 
tions reduce to the Oberbeck-Boussinesq equations (OBE). This approximation, which is 
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formulated on a phenomenological basis, gives quantitatively correct predictions in most 
cases [1]. 

A justification of this approximation, based on introducing a scaling which leaves the 
Rayleigh number invariant is given in [2] (see also [3] and [4]). This approach follows the 
general strategy of taking into account the invariance, under appropriate scaling, of the 
hydro dynamical equations. Such considerations allow, for example, to derive the incom- 
pressible Navier-Stokes equations from the compressible ones, as well as from microscopic 
and kinetic models. 

The main aim of this paper is to derive the OBE starting from the Boltzmann equation 
(BE), which describes gases on the kinetic level, intermediate between the microscopic and 
the macroscopic. To go from the kinetic BE to a hydrodynamical one it is necessary to 
consider situations in which the Knudsen number £, the ratio between the mean free path 
and the size of the slab, is very small. It is well known that the inviscid Euler equations 
correctly describe the behavior of this system for times of order in the limit £ — > [5]. 
To obtain the OBE we need to consider longer times, of order so as to get the effects 
of viscosity and thermal conductivity. To make this possible, we have to study the system 
in the incompressible regime, corresponding to macroscopic velocity fields of order e [6]. 
This scaling would appear, at first sight, to require that the force G be scaled as to be 
consistent with the incompressible regime, i.e. to get a finite force term in the equation 
for the velocity field. However, the case of a conservative force is special from this point 
of view in that larger forces are permissible. In fact we find that G has to be scaled as 
in order to keep the Rayleigh number finite. 

We take the walls to be at fixed temperatures and impose a no-slip boundary condition 
for the velocity field at the hydrodynamic level. This is modeled at the kinetic level 
by assuming that each particle colliding with a wall is reflected with a random velocity, 
distributed according to the equilibrium distribution at the temperature of that wall, i.e. 
we assume Maxwellian boundary conditions. 

Our solution of the BE is given in terms of a truncated expansion in e whose leading 
term in the bulk is a global Maxwellian. The term of order e, denoted by /i, determines 
the hydrodynamic quantities which are close, up to order to the density, velocity 
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and temperature solutions of the OBE. Near the boundary, in a thin layer of size £, the 
hydrodynamical approximation is not correct and we provide a detailed description of the 
solution in this region. The initial datum is chosen to match the expansion up to order e^, 
to avoid a treatment of the initial layer which is more or less standard. 

The validity of the expansion is established up to a time t such that the OBE have a 
sufficiently regular solution by estimating the remainder. This expansion technique goes 
back to Hilbert, Chapmann and Enskog; a rigorous proof of the hydrodynamic limit is 
given for the Euler case in [5] and for the incompressible Navier-Stokes case in [6]. In 
particular, in the absence of gravity, the results of the present paper extend the ones of [6] 
to the case of a fluid confined in a domain with walls at different temperatures modeled 
by Maxwellian boundary conditions. The main technical difficulty for such systems, even 
in the absence of an external force, is in dealing with the terms coming from the boundary 
conditions in the estimate for the remainder. The approach proposed in [5] , and used also 
in [6] , is based on the estimate of some Sobolev norm of the solutions. But in the presence 
of the boundaries this cannot be used because the derivatives of the solution may become 
singular at the boundaries. We avoid the estimates of the derivatives by first looking for 
L2 estimates and then improving them to L^q estimates. This technique was already used 
in [7] and [8] which concern essentially one-dimensional problems, i.e. a compressible fluid 
in a slab with the walls held at fixed temperatures under the action of a force parallel to 
the walls in the stationary regime. Since we consider here the fully three dimensional time 
dependent case we need to modify the method to bound the L^o norm of the remainder 
by using a new method based on a result in [9]. This is presented in Section 4 and in the 
Appendix. A formal expansion including boundary layer corrections was given earlier in 
[10]. 

The case with a force perpendicular to the walls presents extra difficulties stemming from 
the fact that we need good properties of the derivatives with respect to v^, the vertical 
component of the velocity, to get the exponential decay of certain boundary layer terms 
and to control the remainder. Unfortunately, the derivative with respect is singular 
at the boundaries at Vz = 0. To overcome this difficulty we have to decompose the force 
into a part acting in the bulk only and a part acting only near the boundaries, decreasing 
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to zero at a distance of order e. The Milne problem we consider for the boundary layer 
terms involve these short range forces and can be solved by using the result in [11], where 
it is proven that the Vg, derivative is bounded in L2 and L^o norms locally, away from 
the boundaries. This is enough to control the terms appearing in the equation for the 
remainder and in the Milne problems for the higher order corrections. 

While the above results are independent of the nature of the solution of the OBE, we are 
only able to construct stationary solutions to the Boltzmann equation which correspond, 
at the hydrodynamic level, to the behavior of the purely conducting stationary solution of 
OBE. This is due to the fact that the methods in this paper are based on the perturbation 
of a global Maxwellian and therefore require that the temperature difference between the 
plates be small, corresponding to having a small Rayleigh number. We expect to be able 
to construct the purely conductive solution of the Boltzmann equation for any Rayleigh 
number, even when the basic hydrodynamic solution becomes unstable, by perturbing a 
Maxwellian corresponding to the hydrodynamic solution. This involves many technical 
diflBculties and will be discussed in a forthcoming paper. The hope is to extend these 
results to the convective solutions which appear for larger values of the Railegh number. 

2. Hydrodynamic description. 

We consider an incompressible heat-conducting viscous fluid in a horizontal slab A = 
T\ X [—1,1], where T\ is the two dimensional torus of size L. The acceleration of gravity 
is the vector G_ = (0, 0, —G). We specify the temperature on the boundaries z = ±1 as: 

T(l) = T+, T(-l) = T_ = T+ + 6T. (2.1) 

The total mass of the fluid is specified to be m with p = m|A|~^ the corresponding mass 
density. We will assume T_ > T+. Denote 9 = T — T_ the deviation of the temperature 
from T_ and set 

= 0-Ig{1 + z). 
5 

The equations describing the evolution of the velocity and temperature field, u and 9, 
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are the Oberbeck-Boussinesq (OBE) equations (see [12], [4]), which we write as 



divu = 0, 

p{dtu + u ■ Vu) = ijAu — Vp — apOG (2.2) 
^pidtO + u-W9) = kAO. 

Here r\ is the kinematic viscosity coefficient, k the heat conduction coefficient, a = TZ^ 



the coefficient of thermal expansion. 
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p = p —G\\ + z 

and p is the unknown pressure which arises from the incompressibihty constraint. The 
initial conditions are 

u{x,0) =uo{x), Oix,0)^eo{x)-lG{l + z) (2.3) 

5 

for any xE A, with divuo = 0. The boundary conditions for this problem are 

u{x,y,-l,t)^u{x,y,l,t)^0, e{x,y, -l,t) ^ 0, e{x,y,l,t) = -ST - ^G, (2.4) 

o 

for any (x, y) G T'j^ and any positive t. 

Moreover, let r = p — p be the deviation of the density from the homogeneous density 
p, let f = r + pG{l + z)/T_. Then f is determined by the Boussinesq condition 

pV^ + T_Vf = (2.5) 



It can be proved that if uq and are smooth functions of x (i.e. with Sobolev Hs{A)- 
norm finite for some s sufficiently large), then there is t depending on the initial and 
boundary data such that the system (2.2), (2.3), (2.4) has an unique solution, at least as 
smooth as the initial data, for < t < f. We do not give the proof of this statement which 
is rather standard and refer to [4] for details. 

The regime under which the OBE are expected to be valid correspond to a low Mach 
number, a sufficiently weak gravity and a small difference between the temperatures of the 

5 



top and bottom walls. To make this percise we introduce a space scale parameter e and 
rescale the variables as follows: 



(2.6) 



This scaling is a natural one in order to derive the OBE because, under it, the Rayleigh 
number Ra defined as 



is kept fixed. We refer to [2] for the detailed (formal) derivation of (2.2), (2.3), (2.4), (2.5) 
from the compressible Navier-Stokes in the limit e ^ 0, while in next sections we will 
provide its rigorous derivation from the Boltzmann equation. 

We note that our equations do not coincide exactly with the usual Oberbeck-Boussinesq 
(OBE) equations as given in [12], [4] because of the term proportional to G in the boundary 
conditions for 6 and of the quadratic term in G in the definition of p. In the usual 
experimental conditions (see [1]) such terms are much smaller than the others, so one can 
neglect the effect of the variation of the density due to the gravitational force. If we denote 
by ps and Tg the solution of the stationary problem 



with Pg = PsTg and boundary conditions (2.1) the approximation corresponds to setting 
Pg ~ P{1) = p(l)T(l). In this way we would recover the usual OBE. We finally remark 

that the Boussinesq condition (2.5), which is assumed as a "equation of state" in the 
usual discussions of the Boussinesq approximation (see [12], [3]), in our approach is just a 
consequence of the scaling limit (2.6). 



We consider the BE for a gas between parallel planes. We keep the notations of Sect. 2. 
To model the hydrodynamic boundary conditions we choose the so called Maxwellian 
boundary conditions: when a particle hits the walls of the slab (;2 = — 1 or ^ = 1) it is 




^Ps = -Gp. 
az 



AT, = 0, 



3. Kinetic description 
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diffusely reflected with a velocity distributed according to a Maxwellian with zero mean 
velocity and prescribed temperatures T_ and respectively. In the language of kinetic 
theory this means that the accommodation coefficient equals one. The above prescription 
implies the impermeability of the walls, namely no particle flux across the boundary is 
allowed. We introduce as scale parameter e the Knudsen number. The height of the slab 
is 2s~^, hence in rescaled variables z G [—1, 1]. We take for simplicity periodic conditions 
in the x, y direction, and call 

A = {x : (a;,y) e T|,2 e (-1,1)}, VL = {{x,v)\xe K,v eW^} (3.1). 

The BE rescaled according to (2.6) is 

dtf + e-\ -Vr + G- V.r = e-^Q{f, f), 

with 

Q{f,f){x,v,t)= / dv^ / dujB{u>,\v -v^\){f{x,v',t)f{x,v'^,t) - f{x,v,t)f{x,v^,t)} 
Jm^ Js2 

where /S'2 = {a; e |a;^ = 1}, 5 is the differential cross section and v',v'^ are the incoming 
velocities of a collision with outgoing velocities v, and impact parameter uj. We conflne 
ourselves to the collision cross section B{uj,V) = \V ■ uj\ corresponding to hard spheres 
[13]. 

The initial condition is 

r{x,v;0) = f^{x,v), xeA, (3.3). 

The precise form of /q will be specifled below where it will be seen that it cannot be given 
arbitrarily if one wants to avoid a detailed analysis of the initial layer. However, we assume 
the initial datum /o non negative and normalized to the total mass which we set to 1. 
The boundary conditions are: 

f{x,y,-l,v;t) = a-M-{v), v, > 0, t > 0, 

_ (3.4) 

f%x,y,l,v;t)^a+M+{v), v,<0,t>0, 
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(3.2) 



with 

M±(^) = 2^e-V2T,^ (3.5) 

normalized so that \vy\M±{v)dv = 1. The temperature is assumed to satisfy 

T+ =T_(1-2£A) 

with A independent of e, according to the scaUng (2.6). 

The quantities a± must be chosen in such a way that the impermeabihty condition of 
the walls is assured, i.e. 

{v,r) = [ vj'dv = for ^ = ±1, (3.6) 

where we have introduced the notation (/) = J-^s f{v)dv. Condition (3.6) and the normal- 
ization of M± imply: 

a± = ± vj^{x,y,±l,v;t)dv (3.7) 

Namely, a± represent the outgoing (from the fluid) fluxes of mass in the direction z. 
The impermeability condition implies that the normalization of the solution to (3.2) stays 
constant and therefore we will look for solutions to (3.2) which are normalized to 1 as the 
initial datum. 

The macroscopic behavior should be recovered in the limit e going to zero. More 
precisely we expect that for s small the behavior of the solution (3.2) is very close to 
the hydrodynamical one, in the sense that it can be described by a local Maxwellian with 
parameters which differ from constants by terms of order e, and that these terms are 
solution of the OBE. At higher order in e there will both be kinetic and boundary layer 
corrections. Therefore we look for a solution of the form 

7 

r = M + sh + J2 ^"fn + e^R (3.8) 
where M is the global Maxwellian 
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If we put (3.8) in the BE (3.2) we see immediately that /i has to satisfy 

Ch :=2Q(M,/i) = 0, (3.9) 

where C is the linearized Boltzmann operator. (3.9) implies that /i has to be in Null £, 
which means that it is a combination of the collision invariants Mxi with Xi i'^) = ^I'^ii {v^ — 
3T_)/2, for z = 0, i = 1,2,3 and i = A respectively, suitably normalized to form an 
orthonormal set, in L2{M {v)~^ dv) . Hence we have 

4 2 

h = Mj2xMt,x) ^ + ^ + (3.10) 

1=0 P - 

The functions ti{t,x) and/or r,u,9 will satisfy equations to be determined. To write the 
conditions for we decompose them into two parts Bn and 6^, representing the bulk 
and boundary layer corrections. The latter are significantly different from only near the 
boundary. The Bn have to satisfy for n = 2, • • • , 7 

dtBn-2 + v-VBn-i + G-VyBn-2 = 2Q{M,Bn)+ Yl QiBi,Bj) (3.11) 

i+j=n 

where Bq = M and Bi = fi. 

We note that the condition fi = Bi means that there is no boundary layer correc- 
tion to the first order in e. To make this compatible with (3.4) wc need to assume that 
u{x, y, -1, t) = u{x, y, 1, t) = 0, 9{x, y, -1, t) = 0, 9{x, y, l,t) = 2AT_ for any {x, y) e 
and any t > 0. We remark that M + efi, when evaluated for 2; = 1 is not proportional to 
the Maxwellian M+, even with previous assumptions, but differs from it for terms of order 

which will appear in the corrections of higer order. 

To construct the boundary layer terms wc decompose the constant gravity force G_ = 
(0, 0, —G) into three parts: a bulk part Gq and two boundary parts which are different 
from zero in the bulk and near the walls respectively. Their definition is 

G = G+ + G° + G" 

with Gg and G^ smooth functions such that 

r+(.\~! ^' ^-Se<z<l / G, -l<z<-l+5e , . 

0, -l<z<l-25e' 0, -l+26e<z<l ^^'^^^ 
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f: -::?^r (3.13) 

Moreover we have to scale back to microscopic coordinates around z = ±1. Setting = 
£-1(1 z) so that z^ e [0,2£-i] we have that G^{z^) is zero for z^ G [25,2e-\ The 
boundary layer corrections relative to the wall z = ±1, 6^, are chosen to satisfy, for 
n = 2 . . . 5, the equations 



where we have put 

6± = 6± = 0, v = [y^, Vy), V = {dx, dy), = 1, x" = 

£± = 2g(M±,-), AM = £-^[M-M+], M+ = M{p+,0,T+;v), 

and = p + £:r(l). 

Finally the equation for the remainder is 



dtR+-v-VR + G- VvR = \CR + -C^^^R + C^'^'^R + s^QiR, R) + e^A (3.15) 

with 

7 

£Wi2 = 2g(/i,i?), £(2)^^2g(J]£--Vn,i?) (3.16) 

n=2 

and given by 

^ = - ^^(/e + £/7) - • VS7 - i) • - i) • V6f - G • V^(S6 + £^7) 
- (G° + G-) ■ V„[(6+ + £6+)] - (G° + G+) • V„[(66 + £6f )]+ 
2Q(AM,6-)+ J] s'^+^-'QifkJm) 



fc,m>l 
fe+m>8 



(3.17) 



The boundary conditions for these equations have to be chosen in such a way as to 
satisfy (3.4)-(3.6) for f^. Since we are interested in the case = T_(l — 2£A) it is 
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easy to satisfy (3.4) up to the first order in e, because M is already a Maxwellian whose 
temperature and velocity field are chosen to fit with M_, while M + efi is close to be 
proportional to M+ at 2; = 1, up to terms of order s^. 

From the second order on we have to use boundary layer terms to fit boundary condi- 
tions. In fact, as we will see later, the Bn, for n > 2, do not reduce to a^M±. The idea 
is to introduce at one of the boundaries, say z = 1, the correction 6^ so that B2 + 6 J is 
proportional to M_|_ for < 0. On the other hand, the same has to be done at 2 = — 1 
and /2 is modified by 6^ also. This changes again /2 at 2; = 1 by non Maxwellian terms. 
However, since decays exponentially fast, the modification is exponentially small me~^. 
Therefore we impose on the fn the following boundary conditions: 

fnix,v;t) = a~M_{v) (v), z = -1, > 0, t > 

_ (3.18) 

fn{x,v;t) = a+M+{v)+-f+Jv), z = 1, v, < 0, t > 

with the functions 7^^e(i') exponentially small in and such that {'Jn,s'^z) — 0; be 
specified later. Moreover 

a^ = ± Vzfnix,y,±l,v;t)dv (3.19) 

Finally, to fulfil (3.4) we impose the following conditions on R: 

7 

R{x, v; t) = a^M- (v) - ^ £''"^7n,£, z = -1, > 0, t > (3.20) 

7 

R{x,v;t)^a+M+{v)-J2^"~''ln,s, = 1, < 0, t > (3.21) 

n=2 

The initial conditions for -R(x, v; 0) are chosen to be R{x, v;0) = 0, z ^ ±1 for sim- 
plicity. The initial values for the /„'s are partly determined by the procedure below, so 
that only their hydrodynamical part can be assigned arbitrarily. To remove such restric- 
tions one would have to include an analysis of the initial layer, which we skip to make the 
presentation simpler. Finally we impose the conditions 

/ dxdvfn = = / dxdvR (3.22) 
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to ensure the normalization of the solution. Note that this condition on R is satisfied 
because it is true at time t = 0. 

Outline of Solution. 

The equations for the are coupled in a complicated way and have to be solved in 
the proper sequence, which we now outline. The hydro dynamical part of the bulk terms 
is determined by the solvability conditions for (3.11), that we get by multiplying (3.11) 
by Xi)^ = - • - 4, integrating over velocity and using the fact that {Q{f,g)Xi) = 0- The 
solvability condition for (3.11) with n = 2 is 

{Xi[v • V/i + G • V,M]) = 0, z = 0, • • • , 4 (3.23) 

because the Maxwellian M does not depend on x and t. This is equivalent to 

divu = 0, pVd + T_Vr = pG. (3.24) 

The first one is the usual incompressibility condition while the second one becomes the 
Boussinesq condition (2.5), when one defines r = r + p((j'/T_)(l + z). Once (3.24) are 
satisfied, we can deduce from (3.11) with n — 2 the following expression for i?2, where 
denotes the inverse of the restriction of C to the orthogonal of its null space 

4 

+ MY,X^tf\t,x) (3.25) 

The solvability condition for (3.11) with n = 3 is 

{Xi[dth + G ■ V„/i + V ■ VS2]) = 0, z = 0, • • • , 4 (3.26) 

and this produces the equations for u and 0. Let us fix i = 1, 2, 3 in (3.26). Then the first 
term gives the time derivative of pu. The second one reduces to —Gf after integrating by 
parts. Finally we write 

2 2 

{v ® VB2) = {[v®v- \l]B2) + {\lB2) 
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The first term, as is well known, gives rise to the dissipative and transport terms in the 
second of (2.2), while the second one is the second order correction to the pressure P2. 
The result is 

p{dtu + u ■ Vtt) = vAu - VP2 + Gr. 

Using the Boussinesq condition, the definitions of r and 9 and the relation between P2 and 
p of Section 2, we find (2.2)2 ^ Section 2, with r] given by 

2 2 

To get the equation for the temperature, it is convenient to replace X4 in (3.26) by 
X4 = If'i'^ — 5T_)]. A simple computation, using the Boussinesq condition, yields: 

{\[v^ - 5T-)]h) = lp[e - (^G)] + const., 

{v^[v^ - 5T_)]B2) = -KVe + ^-pud{l - z)). 
Collecting the above results we get (2.2) with n given by 

K = {vhv^ - 5T-)C-^[Mvl{v^ - 5T_)]). (3.27) 

Finally, equation (3.26) with i = gives 

a,r = divt(2), t(2) = (4'\tf ,tf). (3.28) 

Summarizing our results so far: we have shown that, assuming u,p,9 satisfy the OBE 
(2.2), (2.3), (2.4), (2.5) up to a time t, the coefiicients U entering in the definition of /i are 
determined. Therefore, once initial and boundary conditions for the OBE are specified, fi 
is completely determined function of {t,x,v). 

On the other hand the hydrodynamic part of B2 is not yet determined, but, by (3.28), 
div t^'^^ is determined in terms of r. Moreover, a combination of tg^^ and t^^^ contributes 
to the pressure p which is determined by the OBE, so that these parameters are not 
independent. 
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The non-hydrodynamic part of B2 depends on the derivatives of r, tt, 6 which are in gen- 
eral different from zero on the boundaries. Therefore B2 violates the boundary conditions 
and we need to introduce to adjust the boundary conditions. We choose 6^ by solving, 
for any t > 0, the Milne problem 

v,—h-e^G-—h = C-h, {v,h) = 

with boundary condition (at z~ = 0) prescribing the incoming flux as the opposite of the 
non hydrodynamic part of B2 at z = —1. The results in [11] tell us that the solution 
approaches, as 2;" —> 00 a function in Null C~ . Thus we set h2 = h — ■, which will 
go to zero at infinity exponentially in z~ and will be the boundary layer correction we are 
looking for. 

In conclusion, we have 
4 

/2(x,^;;0 = M^tf (x,2/,-l;t)x^(^) + 6j(x,2/,2£-i;t)-g2-, ^ = -1, t;, > 0, t>0 

1=0 

(2) 

Since the coefficients tl can be chosen arbitrarily on the boundaries we use them to 
compensate ^2 • To satisfy the impermeability conditions we have to choose = on the 
boundaries. The coefficients of the hydrodynamic part of B2 will, for z 7^ 0, be determined 
by the compatibility condition for n = 4. These are time-dependent non-homogeneous 
Stokes equations (linear second order differential equations) in a slab, together with the 

(2) — (2) 

b. c. — ^215 * = 1) 2,4. Then ' is found up to a constant that is chosen so that the 
total mass associated to /2 vanishes. Finally we get 

/2(x, y, ±1, V, ^0;t) = atM± + 73^,, a± = (x, y, ±1) - gj^^ (0) 

Iterating this procedure it is possible to find all fn- To prove that the terms in the 
expansion have the right properties we use the results in [11] for the solutions of the Milne 
problem with a force, that we state below 

Let F{z) = —Wlz) be a force vanishing at infinity such that V{x) and its derivative 
are bounded. Define 

M = e-^^^)M, L/ = e-^(^)-^2g(VM/,M). 
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Consider the following Milne problem: 



df df J r , I ^ 



< 2; < +00 



lim f{z,Vz) = I < +00 

z— >+oo 



dvv^VMf = 



dvVMs = 



(3.29) 

(3.30) 
(3.31) 

(3.32) 
(3.33) 



In [11] it is proven the following 
Theorem 3.1 

1 ) Suppose that for r > 3 and some u' > there are finite constants ci and ci such that 



sup(l + \v\y\h{v)\ < ci 



sup e'^ ^ sup (1 + |i'|)'^|s(2;, v)\ < C2 



(3.34) 



Then there is a unique solution f e Loo(M+ x M^) to the Milne problem (3.29)-(3.32). 
Moreover there exist constants c and c' such that f verifies the conditions: 



foo e NullL 

sup SM^{l + \v\Y\{f{z,v)-foo{v))\<c 

zeK+ vi 



for any a < c' . 

2) Suppose that for fixed r > 3, £ > 1 and some P > 



sup (1 + 



0^ h ft , I ixT- s 

+ sup sup (1 + |^;|) 



dv7 



(3.35) 
(3.36) 



(3.37) 



for some constant C£ and i ^ 3. Then there are finite constants c and c[ such that 

-d'f d'f. 



sup e^^ sup (1 + It'l)'' [ 



' dvj dvf 



< Cp 



(3.38) 
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for any a < c. 

3) If^ := sup_j,g(o,+oo)[l-^'l + l-^|] exists and is finite then for any S > and for ^ sufficiently 
small there exists a finite constant Cs such that 



sup sup (l + |v|)' 

vei?<* a€(<5,+oo) 



dz 



<Cs, sup sup + 

veR"^ ze[S,+oo) 



df df. 



■ dvz dv. 



<C5 (3.39) 



4. Results in the time-dependent case. 

The main properties of the f^s are summarized in Theorem 4.1 below. 
Theorem 4.1 

Suppose that there is t > such that p{t), ${t) and u{t) are smooth solutions of OBE, with 
\\Vu{t)\\H^ + ||V^(t)||if^ < q for sufficiently large s and < t < t. Then it is possible 
to determine functions f^, n = 2, . . . ,7 satisfying, for < t < t, equation (3.11) and the 
conditions 

fn{x,v-Q)^fl, fn{x,y,±l,v^,Vy,v,^0;t)^a^M±, t > 0, (4.1) 

{A) = (4.2) 
/ dvdxdydzfn = (4.3) 

i]R3xT2x[-l,l] 

Moreover, for any £ > 3 there is a constant c such that: 

sup_|/nk/i < eg (4.4) 
o<t<t 

sup \A\£^h < cq, (4.5) 
o<t<t 

forh< 1/(4T_). Here 

IfU^h^ sup sup(l + |^|)^exp[/i^2]|_^^^^^.^)| ^4_g) 
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Proof 

The proof is achieved by showing that every step of the procedure described in the 
previous section is correct, namely that the conditions on the source and on the boundary 
conditions for the Milne problems that we have to solve at each step are satisfied. Moreover 
we need to check the solvability conditions for the Stokes equations. 

Step 1 

The first step is finding the boundary layer term 63" solving for any (x, y) e the 



Milne problem for 92 — ^2/^ : 

\/^+92{x,y,0,v;t) = B2{x,y,l,v;t)-q^{x,y,v;t), > 0,t > 0, (v^g) = 0, (4.8) 

where M+ = exp[— F+(2;"'')]M+, —dg+V'^ = e'^G'^, B2 is the non-hydrodynamical part 
of B2 given by B2 = L'^ v -Wfi + G- WyM - Q(/i, /i)] . Finally q^{v; t) is the limit 
at infinity of the solution 6^ of the same Milne problem with boundary condition S2, as 
explained in the previous section. 

The force G'^ has been chosen smooth and vanishing as z'^ goes to +00 in such a way 
as to satisfy the assumptions on the force in Theorem 3.1. Furthermore, for s small the 
force term (and its derivative) in (4.7) is small. The boundary conditions verify (3.34) by 
the property of jC~^ (see [13]). Hence, by Theorem 3.1, M^^^^6j satisfy (3.35)-(3.39). 

In the same way we construct 63 imposing the boundary condition in z~ = given by 
B2{x,y,-l,v;t). 

Step 2 

(2) 

As explained above the coefficients t- , z = 1, 2, 4 of the hydro dynamical part of B2 are 
determined by the compatibility condition for n = 4 

{Xi[dtB2 + v-VB:i + G-VyB2]) 

where 



B. = £- 



dth + V-WB2 + G- V„/i - 2g(/i, B2)] +MY,X^ tf^ {x, t) (4.9) 



i=0 
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Proceeding as in the determination of the Boussinesq equation, we find now a set of three 

(2) 

hnear time-dependent non-homogeneous Stokes equations for t] . The non-homogeneous 
terms depend on the third order spatial derivatives of /i. We note that the non linear 
terms in the hydrodynamic equations come from the quadratic term in (3.25), 

while in (4.9) appears a term linear in B2. General theorems for the Stokes equation 
assures the existence of a solution for the chosen boundary and initial conditions. 

Step 3 

Once B2 is completely determined, (4.9) gives the non-hydrodynamical part of B3, B3. 

As before, we introduce the terms bf to compensate S3 on the boundaries z — ±1. The 



term bt is found as a solution of the Milne problem for : J M^g^ — bt 



+ £^G+d,^gs = C+gs + s{x, y, z+, v; t) (4.10) 



with source 



^/^s{x, y, z+, v;t)=v- V6+ + 2Q(AM, b+) + 2Q(/i, b+) (4.11) 
and with boundary condition 



M+93{x, y, 0, v; t) = Bz{x, y, 1, v; t) - {v; t) v^>0,t> 0. 

We have to check that the source satisfies the conditions of Theorem 3.1. 

The condition (3.33) is true due to the properties of Q and to (3.32) for The terms 
of the form Q{f,g) are bounded by means of the Grad estimates 

\M--^Q{f,g)\r-i < C\\M--^f\r\\M-ig\r (4.12) 

so that 

\m:~' Q{fi,b+)\r-i < c\M:yi\r\M:h+\r 

The second term in (4.11) is bounded in the same way 

M-M_ 



\M_'Q{AM, 6+)|,_i< 



M_ 
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\M_^b2\r 



Since T_ =T+(1 + 2£A) we have that k"^^"^^lr-i < + 1) VM|^_i. 

The functions /i and 6^ have bounded norm, hence the third term in (4.11) is bounded. 
To conclude the step we need to bound the first term. It is bounded by using the properties 
for the derivatives with respect to x and y of b2 assured by Theorem 3.1. 

Step 4 

Bs is constructed as in Step 2. Instead the Milne problem for has to be discussed 
since in the source appear also derivatives of 6 J with respect t and v^- We have for 
= &4 that 

^^4 , .2^+ 9 



V 



' " dvz 



^+sa{x, y, z+, v; t) = dtb+ - (G^ + G" )-^b+ +v-Vb+ + 2Q{AM, 6+)+ 

(4.13) 



J2 [2Q{B,,bf)+Q{bf,bf)+Q{bf,bJ) 



i+j=4 

The first term in r.h.s. of (4.13) is bounded by observing that the time derivative of 6^ 
is a solution of the Milne problem we get by differentiating (4.7) with respect to time with 

boundary condition 9(62 (0? v;t), > 0, t > 0. 

The second term in r.h.s. of (4.13) is zero by construction for < z~ < S so that 

dg2 dg2^,\\ 



(G° + G-)V^,6+|^_i < C|l+ sup sup (l + l^^l)' 



and it is bounded because of (3.39), Theorem 3.1, that applies as shown in Step 1. 

Finally the condition (3.33) is satisfied since the velocity flux of the derivatives of 6^ 
with respect to time and velocity is zero. 

The terms /„ with higher n are constructed and bounded in the same way. As a 
consequence, it is easy to see by using the preceding arguments that the term A in (3.15) 
satisfies the bound (4.5) and (4.2). This concludes the proof of Theorem 4.1. 



To complete the construction of a solution to the BE, we have to show that the remainder 
is bounded in norm | • |^^^ defined in (4.6). The remainder has to satisfy (3.15) and the 
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conditions (3.20) and (3.21) with {'~in,e'^z) — 0- Moreover R has to satisfy 

R{x,v; 0) = 0., {v^R) = in ^ = ±1 (4.14) 

which impUes ct^ = ± t!2i?(±l, w; t)(if . To construct the solution of (3. 15), (3. 20), 
(3.21) and (4.14) we first deal with the following linear initial boundary value problem: 
given i:) on X [-1, 1] X M^^ find R such that 

dtR + e~^v -VR + G- VyR = e'^CR + e'^C^R + C'^R + D, (4.15) 

with the same initial and boundary conditions as before. 

Once one obtains good estimates for the solution of this linear problem, the non linear 
problem is solved by simple Banach fixed point arguments, for small e. This allows to 
conclude the existence of the solution and its convergence to the solution of the OBE. 

Solution of Linear problem. 

We consider the linear problem (4.15) with a given D satisfying (vzD) = 0. Put R = 
Therefore the equation for $ is 

dt^ + e-^v • V$ + G • = £-^L$ + £-^L^$ + P$ + (4.16) 

where D = M'^/^D, M = Mexp[£G(^ + 1)/T_] and 

L^f = M-^/^£^M^/^f, j = 1,2 

The boundary and initial conditions are 

^{x,v;0) = 0; ^{x,y,±l,v;t) = a%M±{v)M-^/^ + , v^^O, t>0 (4.17) 

where 

7 » 
C±:=-M-^/2 Ve^-^Tn.e, 4 = ± / v,V^^x,y,±l,v;t)dv (4.18) 

n=2 -^^^^O 

Introduce the L2 norm as 

1/2 



/ 



{ / dxdv\f{x,v;t)\^y 



We now give first the L2 bound for $ and then we provide the L^q bound. 
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Theorem 4.2 



The solution of the linear problem (4-16), satisfy the bound 

Pll < C sup ||L>|| 
teio,t\ 



Proof 

Multiplying (4.16) by $ and integrating on x we have 

+ dxdy^{v,^^){x,y,-l;t)- {v,^'^){x,y,l;t) 



l-dt II $ f = [ dx [ dv 



(4.19) 



To bound the boundary terms in the second line of (4.19), we proceed as follows. We 
consider first the more difficult term 

{v,^^){x,y,l;t)^-{a+f [ dv\v,\MlM-^ + [ dvv,^x,y,l,v;tf. (4.20) 

By the Schwartz inequality, 

= ( / dv\v^\^\v^\^ M^^{x,y,l,v;t)) 

(4.21) 

< / dv\vz\^'^{x,y,l,v\t) I dv\vz\M 

Using the relation between M in z = +1 and M+, the normalization of and the 
relation T+ = T_(l — 2eA), we get, after a straightforward computation. 



dv\vg\M 



2 /r_\i/2 

dvv^M,M~^ = — 



\27tJ 



dvv^M^M'^ 



1 - 2eXy 



(4.22) 



\27tJ (27rT2)2 \2T--T+J (1 - 2£A)2(1 + 2£A)2 

= (1-£2a2)-2<1 + C'£2 

Using this bound and (4.21) in (4.20) we find 



--{v,^\x,y,l,v;t))<Ce [ dv\v,\^\x,y,l,v;t) (4.23). 
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The same argument shows that {vz^'^{x,y, —l,v;t)) is non positive because 



dv\vz\M / dvVzM_M~^ = 1. 

V^<0 JV;,>0 



In Appendix it is proved the foUowing bound for the r.h.s. of (4.23): let ti be any time 
in (O.t]. Then 



e 



(4.24) 



/ dt dxdy / dv\vz\^'^{x,y,l,v;t) 
Jo Jt^ Jv^>0 

< / dt\m-,t)f+e^ sup \\D{-,t)f+ sup |C(-,0I' 

Jo 0<t<ti 0<t<ti 

where |C| = |C+| + IC"!- 

We recall the following crucial properties of the linearized Boltzmann operator L (see 
for example [13]): 

1) 

L = -u + K, 

with K an integral operator and u a smooth function. Moreover, for hard spheres, there 
are two constants vq and vi such that 

+ < ^ix,v) < + \v\); 

2) There is a constant C > such that 

($L$) < -C{u^^). (4.25) 

with the usual decomposition $ = $ + $, with $ and $ the non-hydrodynamical and 
the hydrodynamical part of $ respectively. 

To estimate the operators and wc will use the following estimate on the collision 
operator Q (see for example [14]): for any Maxwellian M and for any ye [—1,1] 

vM J^3 Jk3 
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This inequality and the bounds on the /n's imply the following bounds: 

J dxdv^L^^ < C II II II $ nil M-^/^fi II, (4.27) 

f dxdv^L^^ < C II II II $ nil ^fnW- (4.28) 

n=2 

Note that the presence of the product || ^Jv^ \\\\ $ || depends on from the fact that L} and 
L? are both orthogonal to the collision invariants. 

We integrate (4.19) on time between and ti, and recall that $(a;,f , 0) = 0. With the 
notation $^( •) = $(•, t), we get 



II ^t, f<C / dtl II v^§t 



+ II ^/D^t II II 



e-' II M-VVill + ll M-'/'J2fn 



n=2 



+ II A II' 



(4.29) 



7" 

^0 



+ C / dt||$t|p+£4 sup ||D(-,t)||^+ sup |C( -,*)!' 



0<t<ti 



0<t<ti 



The last line in (4.29) derives from the bound (4.24). The first term in the second line 
is due to the bounds (4.27) and (4.28). Moreover || M~^/'^fi \\< C by the regularity of 
the solutions of the macroscopic equations for < t < t and || M~'^/'^Y^n=2fn ||< by 
Theorem 4.1. 

The following elementary inequality 

-\x^ + (ci£-^ + C2)xy < (ci£ + C2)^yV4 

is valid for any positive e, x, y. We apply it with x = ||v^i>||, y — ||$|| and suitable 
constants ci and C2. We get 



ll'< 

Jo 



dtC 



$(-,t) Ip + ll D{.,t) 



+ 8" sup II f)(.,t) ||2+ sup \Ci-,t)\' 

0<t<ti 0<t<ti 
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In conclusion, by the use of the GronwaU lemma, for e sufficiently small, we get: 



sup II $(• ,t) ||< sup \\D{-,t)\\ + sup |C(-,^)| (4.30) 

Q<t<t t£{0,t\ t£{0,t\ 



Loo hound. 

Let us give some notation: 



dk- 



[x= {x,y,z) : {x,y) e T^,z = ±l} 



S := ^{t,x,v) : t e (0,t],(^,-f;) e 1^} (4.31) 

dS:= |(t,^,v) :t = i),{x,v) e fi} U |(t,x,'f;) :te (0, t], x e aA±, ^ o}. (4.32) 

A and Q, are defined in (3.1). 

We call (pt{x,v) the characteristics of the equation 

dtf+-v-WJ + G-WJ = 

e 



given by 



and define 



(pt(x,v) = (x+-t+^Gt'^,v + Gt). (4.33) 
e 2s 



t~ = inf{t > : (t, ipt{x, v)) e S); t+ = sup{t" <t<t:{s, <ps{x, v)) e ^ Vs < t} 

(4.34) 

Given smooth functions H, i> on S and h in dS, consider the initial boundary value problem 

dtf +^v-VJ + G-V,J- ^vf = H (4.35) 

/(t, X, v) = hit, X, v), it, X, v) e dS. (4.36) 
24 



The solution of this problem is written as 

f{t,x,v) = h{t~,(pt-_t{x,v))ex.p^ - ds^D((ps_t{x,v))^ 

+ dsH(s,(fs-t{x,v))ex.p 1^- J ds' ^i>((fs'-t{x,v))Y 



We introduce the norm 









£ 


P,<1 


/ dx 


7 Mf\' 








L Jm3 J 





and denote by L^'^ the corresponding space. We define the operator as 

^sf = f{s,(Ps-t{x,v))exTp {- ds'^u((ps_t{x,v))y 

We will also omit the apex t when there is no risk of confusion. 

We assume that i> corresponds to the collision rate for hard spheres, i.e. 



Then Nt satisfies the estimate 



\\Nlf\\p,,<CeM- 



p,<i 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



The following Lemma allows to bound the L^'^ norm of Nsf in terms of the L^'^ norm 
of /. It is a generalization to the case with a constant force of the theorem of Ukai and 
Asano [9]. 



Lemma 4.3 



Let 1 < q,p < +00 and a = ^ — ^. Then 



< (— j exp{ 



(4.42) 
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Proof. 

Consider the following change of variables v — > j/: 



1 1 2 1 ?/ ~ 2i 

y = X -\ — VT + —Gt , or V + -Gr — e= 

- e 2e 2 T 



(4.43) 



where we have put r = s — t. We have under the change (4.43) 



d 



Nif\\,,,<r-y-M-^} 



With one more change of variables 



dx 



dy 



/(r + t,2/,£^ + ^GT) 

— T Z 



we get 



Hence the result. 



y-x I 

X — > = S= h -CrT, 

T 2 



This Lemma will be used also for p = +oo and q > d to control the (oo, q')-norm of N^f 
in terms of (g, oo)-norm of the solution of (4.37). 



We write (4.16), (4.17) in the form (4.35), (4.36) with 

H = e-^K^ + e-^K^^ + K^<^ + D :^ e'^K^ + H', 

L is the linear Boltzmann operator such that 

L:=k-u, 

L':=K'-u\ i = l,2. 
Here the u'^ are defined analogously to u as 



(4.44) 



iy^{x,v) = M-^/'^ J dv'doj{v' -v)-uh\f£l, 
iy^{x, v) = M-V^ f dv'doj{v' -v)-u^fn\ 



M. 



n=2 



Finally we set u = u+eC'^ -\-e^iP' and it is immediate to check that v satisfies the assumption 
(4.40). 
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We have the foUowing 



Theorem 4.4 

Let $ be the solution of the problem (4-35), (4-36), with v, H and h given as before. Then, 
for anyq>d, | - ^ < J, 

sup II $ \\oo,g < Ce-'^l'^ sup II b ||2,2 +C£-'^/2 sup || C ||2,2 
0<t<t 0<t<t 0<t<t 

(4.45) 

+ Ce^ SUp_ II B' ||oo,q +C£^ SUp_ II C ||oo,q • 
0<t<t 0<t<t 



Proof. 

By (4.37) we get 

||^||oo,g<||A^^/i||oo,,+^ ds||Ar*iy'||^^^ + £-2^ ds||iV*^$ 
In the last term we substitute to $ its expression (4.37) so that 

II Nlk^ 



|oo,q 



(4.46) 



|oo,g 



<|| N\KNl.h ||oo,5 + / ds' 



NtkmH' \\ooa+e~^ II NtkN',,K^ 



(4.47) 



The boundary term wiU be discussed separately. To estimate the terms containing K we 
use the following Lemma whose proof is given in [9] 



Lemma 4.5 

Let the spaces of functions f{v) such that 

I d^;|/(^;)|"(l + |^;|)^ < oo 
with 7 e R. Let 1 < s < r < oo and riQ = 1 — - -\- -. Then the operator K maps Lt into 

s V y 



We have by (4.41) and by Lemma 4.5 



* ||oo,Q < sup II Nl_h \\oo,q +C sup „ / (isexp{-^^^^^2-^} 

0<t<t 0<t<t Jt- ^ 



/ Is 



(4.48) 
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Using the bound J^* dsexp[—e ^z/o(t — s)] < Ce'^, the estimate for the second term in the 
r.h.s. of (4.48) foUows. We bound the last one by using Lemmas 4.3 and 4.5 as foUows. 



(t 






e 


{t 


-s) 




e 


it 


-s) 




e 



[s-s'). 

e 

.{t-s)J \{s-s'). 



/ e \ '^/i / e vn , ~ 

^<^{W^) ((J37)) exp{-^l(t -.) + (.-.')]} II A'* Ik, 

(4.49) 

To get the first inequahty we have used (4.41) and the Lemma 4.3 which hold for q > d. 
The second step is based on the Lemma 4.5 to replace the L°° norm on the velocity with 
the norm, by taking r = +oo, s = 2, 7 = 0. Then the Lemma 4.3 allows to exchange 
the exponents for space and velocity introducing a factor (e/ {s — s'))^^ exp[—i'Q{s — s') /e^], 
with P = I — ^- Finally the Lemma 4.5 again, with r = q, s = 2 and 7 = 0, gives the 
bound in terms of the norm in space and velocity. 

To get convergence of the s'- integral we need < /3 < ^, so we have to choose q> d , 
i — i < i. From the time integrations we get a factor £4-2(^4+^ ) ^j^]^ ^ = ^' = 
Combined with the prefactor £^l^+^^') it produces £'^~('^+'*'). But /j, + fx' = d/2 so that the 
factor we gain is e^"*^/^. In conclusion 

SUp_ II $ ||oo,g < C SUp_ II Nt_h ||oo,g +C£^ SUp _ || H' \\oo,q 
0<t<t 0<t<t 0<t<t 

(4.50) 

+ Ce-'^/^ sup II $ II 2,2 
o<t<t 

Now we bound the boundary term. We recall that its meaning is as follows: define 

T_ = inf{T > : (fi-rix, v; t) G dS} 

If T_ = t then h = ^{x, v; 0) = 0; if t_ < t then h = ^{x, y, ±1, v; t), in correspondence of 
and {(p-T-{x, v; t))^ = ±1. The boundary conditions on $ are given by (4.17) with 



a% = ± Vg^/^^{x,y,±l,v;t)dv 
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Equation (4.37) allows to express in terms of y, —1, v, t), Vz < and H. In fact 
by (4.37) 



^{x,y,l,v;t) = h{t , (^t-_t(a;, y, 1, v)) exp | - ^ ds^iy{(ps-t{x,y,l,v))^ 
+ dsH{s,(fs-t{x,y,l,v))exp^^- ds'^v[(fs,_t{x,y,l, 



(4.51) 



We remark that the characteristic which is on the boundary z = 1 with velocity 
V such that f ^ > at time t~ — t had to start either from some point in the bulk or 
from the boundary z = —1. In the first case we have /i = in (4.51) otherwise we get 
h = ^{x,y, —l,v;t) for Vz > 0. In the latter case it is shown in the Appendix that the 
integral dsiy^ips-tixjy, —l,v)) is bounded from below by Ce. As a consequence, we 
can estimate the exponential in the first row of (4.51) as e~*-^/^. 

We exploit the same argument to deal with and use (4.37) to represent $ on the 
boundary z = —1. There is a difference with respect to the previous case: due to the 
presence of the force the characteristic (p which is on the boundary z = —1 with velocity 
v,Vz < at time t~ — t can start from the bulk, from the boundary z = 1 with negative 
Vz or from the boundary z = —1 with positive Vz- In the Appendix it is shown that in the 
latter case the exponential factor in the first row of (4.51) allows to gain a factor Ce^ . 

We discuss explicitely the bound for a^. The case of is dealt with in the same way. 
By using the representation (4.37) we get 



sup \a^{x,y,—l;t)\= sup / Vz^sfM^{x,y,—l,v;t)dv 
< C sup / dv\<^\'^{x,y,-l,v]t) < 

■Jv, 



(x,j/)eT2 



C sup [e ^\a+\+e^\aji\]+e^ \\C\\q+ I ds\\NlH\\ < 

{x,y)eT^ Jt- 

C£^||$|| +/ ds\\NlH'\\ +e-^ I ds\\Nlk^\\^q 



(4.52) 
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The first bound is obtained from the inequahty 



sup |ct^| < sup 



KeT2 X [ 



3 / 



M\^{x,v;t)\dv 



< C sup 

kGT2x[-1,1] 



dv\^\''{x, v; t) 



1/9 



\oo,q 



As explained before, to get a bound of the last term in (4.52) in terms of the L2 norm we 
need to iterate the procedure and use again the representation formula (4.37) in the last 
term 



ds II Ntk^ 



< e~ 



ds II NlKN^.h lloo.g 



+ 



ds I ds' 



NlKNt,H' \\oo,q II NlkNt,k^ 



\oo,q 



(4.53) 



In this way all the terms we got in (4.53) are analogous to terms already discussed in the 
first part of the proof but the term containing h. The problem with this term is that it 
can be estimated in terms of the Lqo norm of $ but we need to gain a small factor and 
cannot provide it. Hence we have to repeat the previous argument and use (4.37) 
to represent $ on the boundary in terms of the function evaluated on the point on the 
boundary reached after a finite amount of time. 
We get 



IKL + 



/; 

ds II NlKNl. [M-i/2[M+a+ + a^M"]] ||oo,g + 

NlkNl.Nl_H' \\oo,q II NlkNl.Nl_k^ ||oo,g 



e - j ds\\NlKNl.h,,oo,q 



(4.54) 



^ " I ds J dr 



By using the properties of Ng and k and the bounds on the integral in the exponential as 
discussed before the second term in (4.54) is bounded by £^ || $1 



\oo,q 
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Finally we get 



sup \a%\<Ce^\\^\\^^^+e^\\C\U+ I ds\\NsH'\\_ ^+ 



\oo,q 



NlKNI,H' ||oo,Q +s~^ II NtKN',,K^ 



giv^.g/iv ^ ||00,Q 



(4.55) 



Now we can apply all the previous argument to get the following estimate for the boundary 
term in (4.50) 



|A^t-/i ||oo,</< sup [|aj| + |a^|]< 

Ce^ II C ||g +C£'sup II H' ||oo,<? +C£-'^/2sup II $ ||2,2 



(4.56) 



Finally by (4.30) we have 



sup II $ \\oo,q < Ce'^ sup II C llg +C£^ SUp || H' 
0<t<t ' 0<t<t 0<t<t 



1 00,17 



+ Ce-'^/^ sup II D ||2,2 +e-'^/^ sup || C II2 

0<t<t 0<t<t 



(4.57) 



This concludes the proof of Theorem 4.4. 

To get the bound for $ we need to estimate the || • || 00,00 norm in terms of the 
II • 1 1 2,2 of $. The use Lemma 4.5 allows to prove the following 

Theorem 4.6 



Define 



I/I, = sup sup |(1 + |^|)VI- 

xeA vER^ 



Let $ be the solution of (4. 16), (4.17). Then if r > 3 



sup i^ir ^ 
0<t<i 



sup j^l^ + sup IClr 
0<t<i 0<t<i 
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Proof 

By (4.37) we have 



\^\o <\Nt-h\o + ds\NsH'\Q + e-'^ ds\Nsk^o 
By (4.41) and Lemma 4.5 (used with r = +oo and s replaced by q), ior q > ^ we have 



l^lo < C'e^ sup \h\o + C sup \H'\o + C sup \^\oo,q 

0<t<t 0<t<i 0<t<t 



By (4.57) we get 



|$|o < C SWp_\h\o + Ce SUp_|i?'|o+ SUp_ || f/"' ||oo,g +C£ SUp_ || C ||oo,q 
0<t<t V 0<t<t 0<t<t 0<t<t 



+ Ce"^/^ ( sup II D ||2,2 + sup II C II2 

\0<t<t 0<t<t 

Because of the factor (1 + |)'~ in the definition of | • and of the boundedness of the space 
domain, if r > 3 we have 

II / ||oo,,< C|/|„ ||/||2,2<C|/|,. 

Hence 

sup |$|o<C sup \h\o + Ce~'^/'^l sup \D\r + sup |C|r)+C'£^ sup (4.58) 

0<t<t 0<t<t \0<t<t 0<t<t J 0<t<t 

The boundary term |/i|o < C'|ckJ|oo + Ick^Ioo + iCloo has been estimated before by using 
(4.56). The term containing H' will be bounded using the Grad estimate 

\M-'/^Q+{f,g)\r < C\M-'/y\r\M-'/^g\r, (4.59) 

where Q"*" is the gain term of the collision operator. 

By using (4.59) and the definitions of K*, i = 1, 2 we get 

II \\r< C II M-^/^fi \\r\\ $ llr 
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II K^^ \\r<C\\ 11,11 $11,. 

n=2 

Finally the regularity property of the hydrodynamic solution and Theorem 4.1 give 

II ||r< II $ ||„ II ||,< C II $ II, (4.60) 

Hence, for e small 

sup |$|o < C'£-'^/2| 1^1^^ \ (4_gl) 

0<t<t \0<t<t 0<t<t J 

Finally we improve the | • |o norm to | • |, norm by means of the Grad estimate and the 
representation (4.37) of $ 



|$|, < C £^ sup |-H''|r+ sup ICjr + SUp |$|r-l 

y o<t<t o<t<t o<t<t 

By (4.59), iterating the previous inequality we get, for e small, 

sup < Ce^ sup_|L>|, + C sup_|$|o + C sup JCjr- (4.62) 

t 0<t<t 0<t<t 0<t<t 

By putting together (4.61) and (4.62) we eventually get 



sup |$|, <C£~'^/M sup |£>|, + C sup Id 

0<t<t \ 0<t<t 0<t<t 



r 



SO proving the Theorem. 

Non linear case 

We conclude our discussion by proving the following 

Theorem 4.7 

There is Sq such that, if e < Eq, for < t < t there is a unique solution to the initial 
boundary value problem (3.15), (3.20), (3.21) verifying the following: for any positive 
integer £ there is a constant c > such that 

\R\i,h < Ce (4.63) 
33 



for any h < 1/(4T_). 
Proof 

The equation we have to solve is 



with 

and A given by (3.17). 

By Theorem 4.6, for d = 3 we have that 



sup 

0<t<i 



\^\r<e^''^C\ sup \M-^/'^Q{\f^^,\f^^)\r+ sup |M-^/2^| 
yo<t<t o<t<t 



+ sup ICIr 

0<t<t 



By (4.12) 



sup |$|^ < Ce^/M ( sup |$|^)V£sup|M-^/2A| J +C£-2/2g^p|^|^ (^4_g4) 
o<t<t y o<t<t 1/2 y 1/2 

By the arguments in [7] we then get 

sup < S^/^C sup \M-^/^A\r + Ce'^/^ sup ICIr 

t 0<t<t 0<t<t 

The term ( decays exponentially fast in e. Hence, by Theorem 4.1 the estimate (4.63) 
follows. 

5. The stationary case 

The main difference between the results for the time dependent case of previous Section 
and those for the stationary case we are going to present is in the restriction to small values 
of the Rayleigh number, we need to deal with the stationary problem. Therefore, at the 
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hydro dynamical level we are confined to the purely conductive solution. We hope to be 
able to extend our method to the convective solutions which appear for larger values of 
the Rayleigh number. The proof follows by argument quite similar to those presented in 
[7], [8] to which we refer the reader for more details. In this section we only give a sketch 
of the proof. 

We start by recalling the stationary setup. We look for one-dimensional solutions, 
namely for solutions not depending on x and y so that the equation we consider is 

v,d,r + eG ■ V„r = e-^Q{r, f) (5.1) 

The boundary conditions are: 

r{-l,v) = aJM_{v), V, >0, (5.2) 

r{l,v)^a+M+{v), v,<Q, (5.3) 

with M±{v) given by (3.5) and a± now independent on a;, y and t, given by (3.7), so that 
satisfies (3.6). 

We construct the solution in the form (3.8) with to be determined according to a 
bulk and boundary layer expansion. These terms are computed as in the time-dependent 
case and a theorem similar to Theorem 4.1 can be stated also in this case. We only discuss 
the remainder equation because its solution requires a different technique. This equation 
has the form 

—R - eG—R = -CR + C^^^R + eC^'^'^R + e^QiR, R) + e^A (5.4) 
oz avz e 

with £(1) and C^'^^R defined in (3.16). Moreover, A is given by 

(G'0 + G'+)A[(6-+£6-)]+ e'^+^-'QifkJm). ^^'^^ 

k,m>l 
k+m>8 

The boundary conditions on R are given by (3.20), (3.21). R satisfies the normalization 
condition (3.22) and the vanishing fiiix condition 

dvV;,R{z,v) = for [-1,1]. (5.6) 
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The theorem below summarizes the results about the existence of stationary solutions. 
Theorem 5.1 

Let M be the Maxwellian with parameters p, T_ and vanishing mean velocity. Put 

with f and 9 the thermal conduction solution of the OBE corresponding to the temperatures 
T- and — T-{1 — 2eX), namely 

9 = T^X{l + z), r = -p[^-X]z. 

Then there are Aq > and £o > such that, if X < Xq and s < Eq, there is a stationary 
solution to the boundary value problem above such that 

\f^-(M + efi)\r<Cs^ (5.8) 



Sketch of the proof. 

We follow the strategy of the previous section: first we get an L2 bound and then the 
bound. In the present case we cannot use the initial condition to satisfy the normalization 
condition. Therefore, we satisfy the conditions (3.22) and (5.6) by choosing the constants 
along the lines of [8]. 

Observe that (5.6) is satisfied for any z G [—1,1], once it is satisfied at one point, because 
J dvVzR{z, v) does not depend on z in consequence of (5.4). 
We write R as 

i? = /($)M+V^$ (5.9) 

with 

I{^) = - f dz f dvR{z,v), (5.10) 

J-l Jr3 

SO that (3.22) is satisfied. Therefore we have — (T_/27r)~^/^p~^/($). It is easy to 
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check (see [8]) that the function $ has to solve the foUowing boundary value problem: 

oz ovz e 

$(-l,^) = M-VV V, >0, 

_ (5.11) 
$(1,-?;) =/3flM+(v)M-i/2 + C+M-i/2 v^<0 



dvv^^M^/'^ = 

where = ~ Sn=i ^"^~^7n,e Pr — '^r ~ '^R- '^^^ linear operator J\f^ is defined by 

6 

= £L2$ + /(Vi^$)[^£"L/„ . (5.12) 

The non linear term is given by 

Q($, $) = -^g(VM$, Vm$) + 2/( V^$)L$. (5.13) 



In this way there is no normalization condition on the function The quantity 
represents both the outgoing flux oi fa at z = —1 and the integral of R over z and v. The 
constant is determined so that R satisfies condition (5.6) at the point z = 1, i.e. 

pR = [ dvv,R{l,v)+ [ dvv.C. (5.14) 

To construct the solution of (5.11), we first consider the following linear boundary value 
problem: given D on [— 1, 1] x and on {v e M.^s.t. Vg ^0}, find R satisfying 

1 ~ ~, 

v,— -eG— = -L^ + L^^+M^ + D, (5.15) 

oz OVz s 

and the last three conditions (5.11). 

As usual we introduce $ and $ the hydrodynamic and the non-hydrodynamic part of 
$ respectively. Multiplying (5.15) by $ and integrating on velocities, we have 

^ S 

By integrating over z, using (4.25), (4.27) and (4.28) we get 

l{v,^^){l)-l{v,^''){-l) < -C- II $ IP +C[X+e] II $ II [ II $ II + II $ II 1+ II II II $ II . 
2 2 e 
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Here || $ |p= /^^ dz J dvcfP'v. One can check (see [8]) that the l.h.s. of this inequahty is 
positive. Therefore we get 

II $ ll^< Ce\ II $ II II $ II +£ II II II $ II 



Using the inequahty xy < kx^ + y'^/4:k, with a;=||$||,j/ = £||$|| and a suitably smaU 
k, we find 

II $ f < Ce^A II $ f +£ II II II II +C£2 II i) f (5.16) 



To bound the hydrodynamical part, we multiply (5.15) by v^^i, = V Mxi, ^ 7^ 2 
and integrate over [—1,2;] x R^. 

Denoting Pi{z) —< Vz'^i^ >, we get 



Piiz) =pi{-l)+ / dz' / d'W'y^^' 



1 = . 



L$ + + £A/'$ + £G^$ 

£ 



Following [7], pag. 68-69, we can prove that pi have the following estimate: 



|j3i(-l)| < \CX II $ II II $ II +- II II II $ II 1 

£ 



Using (5.16) this inequality becomes 



1 



Pi{z) = [CX II $ II II $ II +- II L> II II $ II ] ^ + C 



- II $ II +(A + £) II $ 

e 



(5.17) 



Since {vz^V M) = 0, we can decompose $ as 



Therefore 



with B a non-singular matrix. This allows to estimate hi and hence ^ as 

$ ||< C 

-£ " 
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- II $ II +- II D II +A II $ 

£ £ 



(5.18) 



so that for A small enough we get 

II ^ ll< c 
Combining (5.16) and (5.19) we get 



- 11*11 +i \\D 

e e 



(5.19) 



II $ ||< C II D II, II d ||< C- II D II (5.20) 

Note that the only point where we need A small is to pass from (5.18) to (5.19). 
L°° bounds. 

To write the integral form of the linear equation (5.15) we follow the approach in [11]. 
The notation is as follows: 

The "total energy" is E{z,v) = vl/2 + V{z) with V{z) = eG{z + 1). The lines with fixed 
E are the characteristic curves of the equation 



df df 

OZ OVz 



For E{z,v) > V{z') we define 



vsiv, z, z') = y/2E{z,v)-2V{z') 

v{z, z') = {Va:,Vy,Vs{v,Z,z')) 

TZv = {vx, Vy, -Vz); TZv{z, z') = {vx, Vy, -V3{v, x, x')) 
Moreover call z~^ {z, v) the function implicitly defined by the equation 

E(z,v) = V{z+) 

Finally put 

Jz' V3{v,z,r]) J^, V3{v,z,r]) 

Consider the equation 

a$ a$ 1 1 

v,^-eG^ + -D^^-H (5.21) 

OZ OVz £ £ 
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with boundary conditions 

f{±l,v) = h{v)^, v,^0 
The solution of (5.21) can be written in an integral form as: 
for Vz > 0: 



f{z,v) = h {v{z,-l))exp~nz-i{v) 

+ / dz' -H{z',v{z,z'))exp — ^z,z'{v), 



1 ev3{v,z,z')' 



for Vz <0 and E < V{1): 



f{z, v) = h {v{z, -1)) exp -\^Vlz+-i + ^Ti^z+, 



+ 



/ dz'—f 7kH{z', {v{z, z'))) exp / + -7ei]^+, 



+ j dz'-—^ -H{z',Rv{z,z'))exp--nnz,z' 

Jz ev3{v,z,z') e 



for < and £; > V{1): 



1 



f{z,v) =ex.p--ni^z{v)h+{v{z,l)) 



+ 



Jz £ 



— — -izH{z', v{z, z')) exp --Qz',z- 

ev3[v,z,z') s 



We can write the previous formulas in a compact form as 

f = Vh + TH 



where 



Vif{z,v) = x{vz > 0)h-{v{z,-l))exp~nz,-i{v), 

V2fiz,v) = xivz < 0)x{E < Vil))h-{viz,-l))exp-^[Qz+,-i + nnz+ 

Vsfiz, v) = x{vz < 0)X{E > V{l))h+{v{z, 1)) exp -^l^i,, 

3 



40 



The definition of T is given in a similar way. 

By slightly modifying the proof in [11] to take into account the factor e it is possible to 
prove the following Lemmas (see also [7]). 
Lemma 5.2 

For any integer r > there is a constant c such that the integral operator T satisfies the 
following inequality, 

IT H\r < c 



H 

V 



(5.26) 



Here 



|/|, = sup(l + |^|)'-|/(^,^) 



Lemma 5.3 

For any 5 > and for any r > 2 there is a constant Cg such that 



N{T H) < -^Cs II ly-^/^H \\ +S\H\r. (5.27) 

Here 

iV(/)= sup (I dv\fiz,v)f) . 
ze[-i,i] \^R3 / 

We write (5.14) in the form (5.21) with 

H = k^ + eL^ + eL^^ + eU^ + eD (5.28) 

and = ^(±1, v) given by (5.11). Combining these Lemmas and using the properties of 
the operator L one gets 

Theorem 5.4 

There exists a constant C such that for any r > 3 the solution of (5.15) verifies 

< II ly-^/^^ II2) + ^ II siy-^/'^D \\ +e\D\r + \Vh\r (5.29) 
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Noting that || v ^I'^D ||< \v ^-Djs and using (5.19) we get 



l^lr < C—j=.\D\r-X + \Ah\, 



(5.30) 



Finally (5.30) implies 



\RM-^I\ < Cy/e\M-^/'^D\r-i + \Vh\, 



(5.31) 



The term containing h involves Pr which still depends on $. To estimate it, we follow 
the method in [8] . 

Equation (5.25) allows to express Pr in terms of TH and the restriction of $(— l,v) to 
Vz > 0. We have the estimate 



/ 

Jv, 



v,My'hil,v)\^ \ I v,M:''{h-{v{z,-l))e^p--n,,-i{v) 

>0 Jv^>0 ^ 

+ I dz'-^^—-H{z',v{z,z'))exp--n,^,,{v),}\ 



rl/2 



1 



Li svs(v,z,z')' 
<\s\ + \ I v,M^/^TH\, 



(5.32) 



with s± = C^M ^/"^ and \s±\ = sup^^^Q |s±(v)|. By (5.25), using the Schwartz inequality 
and 

/■^ 1 1 1 f°° 

/ dz 7 --^(2, f (1, ^)) exp — Qi z(v) < C / dzexp\—z]<C, 

J_l £^3(^,1,^) £ ' Jo 



we get 



Jv,>0 Jv,>0 V3{V,1,Z,) 



1/2 



<e-'/^C \vz\'/'M'/^\ dyu-'H 



Vz>0 



-1 



1/2 



Finally, using again the Schwartz inequality and recalling the expression of /3_r in (5.14) 
we get 



/3R<c[e-'/^ II ly-^H II2 +|s| 



(5.33) 
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Using now the form of H in (5.28) 



(3r < c[£-l/2 II ^-1$ 11^ ^ II 11^ ^|^|]_ 

As a consequence (5.31) becomes 

\RM-^/^\r < C^/e\M-^/^D\r-i + C\M-^/^dr + \M-^/^C\r (5.34) 

which concludes the analysis of the linear case. The non-linear problem is dealt with by a 
fixed point argument and the final result is 

Theorem 5.5 

There exist A > 0, eo > and a constant C such that, if X < Xq and e < Sq, there is a 
solution to the boundary value problem (5.4), (3.20), (3.21), (5.6) such that for any r > 

|M-V2i?|^ < C£3/2|^-l/2^|^ (5.35) 

This concludes the proof of Theorem 5.1. 
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Appendix 

In this appendix we show how to get L2 estimates for the boundary terms. We will 
prove (4.24). To this end it is enough to prove 

S+:=s[ dt [ dxdy [ dvv^\^\^{x,y,l,v;t)<e'^[ dt \\ ^ f + sup C(-,^)l^ 
Jo Jt^ Jv^>o Jo o<t<ti 

(Al) 

where Ht{-) := H{- ,t) is defined in (4.44). In fact by substituting the expression of 
we get 

6^ II i^-'/^Ht f<\\ v-^l^k^t f II v-^/^K^t f 

+ II v-^/^k\t f II Dt f< C \\ \? \\ Dt f 

so proving (4.24). 

We use (4.37) to express the value of the function $ in the point z =1. We remark that 
the characteristic (p which is on the boundary z = 1 with velocity f , f ^ > at time t~ — t 
had to start at time either from some point in the bulk or from the boundary z = —1. 
In the first case the boundary term h in (4.37) is zero otherwise we get $(a;', y', —1, v'; t), 
v'g > where 

x' — X + Vx{t~ —t), y' — X + Vy(t~ — t), v' = V — G{t~ — t). 

Hence 

\^\'^{x, y, 1, v; t) < 2|$|2(t-, x', y', -1, ^;')) exp | - 2 ds^i^{ips-t{x, y, 1, v))^ 

dsH{s,(ps-t{x,y,l,v))exp^^- ds' ^^{(Ps'-tix.y.l.v))^ 

{A.2) 

By the boundary condition on $ we have 

^{t-,x', y', -1, v')) = a^{t-,x', y', -l)M_{v')M-'/\v') + C {v'), < > 
so that by Schwartz inequality we get 

/ v'Jv'\^\^{t-,x',y',-l,v'))<C{a^)\t-,x',y',-l) + C [ v,dv\C\^ 

Jv'^>0 Jv^>0 
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+ 2 



But ctjj is expressed again in term of the value of $ on the boundary z = —1 and for 
negative ^-component of the velocity, namely 



M^t-,x',y',-l,v) 



Vz<0 



< c 



\v^\ ,a;',y',-l,v) 



v,<0 



We observe now that the integral 



dsu{(ps_t{x,y,-l,v)) 



when evaluated on the characteristics going from one boundary to another is bounded 
from below. In fact we can use the bound on u , i/(x,v) > Co(l + \v\) to check that this 
integral for the trajectory perpendicular to the boundary, namely Vx = Vy = 0,Vz > is 
greater than e times the width of the slab and this is the worst case since for the other 
trajectories is even greater. As a consequence, we can estimate the exponential in the first 
row of (A. 2) as e~^l^ . In conclusion we get 



dt J dxdy J dv\vz\\^f{t , x', y', -1, -y')) exp | - 2 ^ ds^iy{(ps-t{x,y,l,v))^ 



/ dt dxdy / dv\vz\ \^\'^ (t, x,y, —l,v)) + sup || ^ 

Jo J JV;,<0 0<t<tl 



(A.S) 



To get the previous expression we have replaced x', y' and t~ with x, y and t, since the 
respective Jacobians are equal to one. 

In this way we still do not have any explicit estimate of the boundary term in (A. 2), 
because it contains the values of $ at 2; = —1 for negative which is still unknown. 
However, we can use the representation (4.37) again to express it back in terms of the 
function in 2; = 1 so to get a set of coupled equations for the boundary terms. Hence we 
consider 



JT'^ 



—el dt j dxdy I dvvz\^\ {x,y, —l,v;t) 



v.<0 



Equation (A. 3) implies 



S+ < Ce-^/'[S- + sup II C" 11^] + £n+ 

0<t<ti 



(AA) 



{A.5) 
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where 



dt I dxdy I dv\vz\ I dsH (^s,(ps-t{x,y,±l,v)) 
Jt^ J Jt- 



expj-^ ds'^u{(Ps'-t{x,y,±l,v))^ 



As above we use (4.37) to find a bound for S . Define 

l^^l^ 1 



E 



2£2 £ 



+ -G{z + 1). 



We have: 



for <0, E>0 

^{x, y, -1, v; t) = ^{t~,x', y', 1, v')) exp | - ^ ds^u{(ps-t{x, y, -1, v)) | 
+ dsH{s,(fs_t{x,y,-l,v))expi^- ds' ^iy{(fis'_tix,y, -l,v))y 

For <0, E <0 

^{x,y,-l,v;t) ^ $(t",a;",|/",-l,t;"))exp| - ds^iy{^'^_t{x,y, -l,v))^ 



(A.7) 



where (x", y", —1, v") — j/, —1, v) and (/p' is the characteristics starting from —1 

with v'l >Q. 

These two formulas correspond to the case in which the characteristics starts from 
z = 1 with vefocity < and to the case in which it starts from z = —1 with v" > 
and then come back to the boundary z = —1. This second possibihty appears because of 
the presence of the force: the kinetic energy is such that the trajectory does not reach the 
opposite boundary but instead goes back after a time depending on the balance between 
kinetic energy and potential energy. In that case the total length £ for a trajectory with 
Vx = = Vy is given by 

^= ^[<f • 
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The other trajectories give a larger value to this integral. 

Hence the square of the second term in the first row of (A. 7), integrated on time, space 
and velocity, is bounded as 



B 



:= f'dt I dxdy [ \v':\^\\x",y",-l,v";t)eM-^vl} 
Jo Jt^ Jv^<o ^ 

f f C 

dt / dxdy / \vz\^\'^{x",y",-l,v'J.,Vy,Vz;t)exp{ — -vj}, 

Jt^ Jv;,>o ^ 



(A8) 



using that v'^ = —Vg. The function $ in the integral in (A. 8) is relative to 2; = —1 and 
> so that we can express it in terms of as follows, for Vz > 0: 



< C 



dv\vz \VM^x", y\ -1, vl, v'' V,; t) 



M + C 



<CM / dv\vzm'{x'\y",-l,v'^,v'^,Vz;t) + {Cy 



{A.9) 



Equation (A. 9) implies for B the following bound 



B<C [ dt ! dxdy I dv\vz\\^'^{x",y",-l,v'l,v'li,Vz;t) 

Jo J Jv^<0 

/C f c 

dv\vz\Mexp{-^v'^J + C / dv\vz II C"Pexp{-— ^;^} 



(^.10) 



We have that 



/ 



C 

dv\vz\Mexp{ — ^v^} < Ce 



(All) 



Finally 



'C r dt I dxdy I dv\v,m\x" ,y" ,-l,v'l,v';,Vz-t) + sup || C 

^0 J Jv^<0 0<t<ti 



B <e 



= CeS~ + sup || C 
t 



- Il2 
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The contribute of the first term in (A. 7) to S is 



f dt f dxdy f dv\vz\\^\'^{x',y',l,v';t ) exp | - /" ds\u{(ps_t{x,y, -l,v))\ < 

Jo J Jv^<0 ^ Jt- ^ 

8 dt [ dxdy I dv\v,\[M\alWx\y\l;t-) + 

Jo J Jv^<0 

xexp|-y^ ds^{ips-t{x,y,-l,v))^ < 

I dt I dxdy I (if |vz||$p(a;',y', + sup || C"^ ||' 

Jo J Jv^>0 0<t<ti 



= Ce-^/'[S++e sup || C+ ||^] 

0<t<ti 



Summarizing we can write 



{A.12) 



S- < Ce[S- + exp[-C/£]5+ + e sup || C" f ] + sH {A.13) 

0<t<ti 



where 7i — 7i~^ + H . This imphes, for e smaU, 



S- < Ce-^/'[S+ + e sup || C f ] + eU {AAA) 

0<t<ti 



Putting together (A. 6) and (A.13) we get 

S+<c{en+ sup II cf) 
V o<t<ti / 

s-<c{en+ sup iicil') 



(A.15) 



To conclude the argument we now provide an estimate for Ti. By the Schwartz inequahty 
for the integral on ds, we have 

1^"^! < j dt j dxdy J dv\vz\ ds\^\^ (s,ips-t{x,y,±l,v)) 

dsv[s,ipa-t{x,y,±l,v))eicp^~ ds'^v{ips'-t{x,y,±l,v))^ < {AAQ) 

<Ce^ j dt j dxdy j dv\vz\ j ds\-j=^ {s,ips-t{x,y,^^,v)) 
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Consider the term 



:= J dt ds J dv\vz\\^\^{s,(pa-t{x,y,±l,v)) 

and the foUowing change of variables {t,Vz) i^jw) 

e = m ±1 + -Vz{t -s) + ^G{t - sf, w = v,+G{t- s) 
e 2e 

whose Jacobian is elf^l"^. Denote also by t{^) the inverse of ^{t). We have 

V ■.= e di I ds dw\^\ {s,x{^,w),y{^,w),^,Vy,Va;,w)) 

Jm Jt-(£) J 



Hence 



tl 

JQ 



\n\<e^ I ds I dx I dv\^\'^{s,x,v). (A17) 



which implies the bound (A.l). 
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